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Abstract — This paper presents an innovative approach for the analysis of the Power Quality 
Disturbances both qualitatively and quantitatively. The proposed method employs the phase 
shifting property of Hilbert Transform for the accurate detection and computation of the 
characteristic magnitudes of the power quality disturbances along with the time of their 
occurrence. This facilitates for the real time detection and characterization of various 
disturbances such as voltage swells, voltage sags, voltage fluctuation, harmonics and 
transient oscillation accurately. The various disturbances have been simulated on the 
LabVIEW platform and the phase shifting property of Hilbert Transform has given 
satisfactory results. Real Time Signals were generated and digitalized by the aid of Data 
Acquisition (DAQ) card, which were processed in the LabVIEW environment, to yield 
immaculate results indicating the characteristic magnitudes and time of occurrence of 
disturbances. 

Index Terms — Power Quality Monitoring, disturbance detection, Hilbert transform, phase 
shifting, LabVIEW 

L Introduction 

The introduction of Power Electronics equipment in power systems and also the increased usage of various 
time varying and typical loads at the consumer end has resulted in disturbances in the Power Quality which 
are diverse and significant. The disturbances broadly fall in the category of voltage swells, voltage sags, 
voltage fluctuation, harmonics, inter-harmonics and transient oscillation. The research undertaken in power 
quality till date was successful in providing tools to detect and monitor each of the disturbances discretely 
and individually, but no effective common technique has been formulated to accurately detect and classify 
the disturbances in real time [l]-[2]. 

The typical detection methods involved are based on time domain methods or frequency domain methods 
making them suitable to detect only a characteristic type of disturbance at a given instant. The proposed 
method of application of phase shifting property of Hilbert transform is inherently a steady inter-convertible 
from time and frequency domain making it an efficient technique for the detection of disturbances of a wide 
range with uncharacteristic properties. The method even provides with the time of occurrence of the 
disturbance making it suitable for real time application. Coupling the same with LabVIEW platform and 
executing the method has given satisfactory results. 
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II. HlLBERT TRANSFROM OVERVIEW 

1 

The HT of the function x(t) is defined as the convolution between x(t) and h(t) = — . The expression is 
H[x (t)] = x h (t) = x(t)*h(t)= x(t)* — 
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~ 1 ^ (1) 
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In formula (1), * represent convolution; H[x(t)] and x h (t) represent a HT in time-domain. Different from other 
transform methods, the HT of time-domain function is still in time domain, the HT in frequency domain can 
be expressed as 

X k (f) = X(f)H(f) (2) 

Because the FT of h(t) , 

F[h(t)} = H(f) = -j sgn(/) = ~/{ \ f > ° (Sgn is the sign function ), 

I, / 

So as long as multiply the negative frequency of the X(f) by j and multiply positive frequency of the X(f) by 
-j , it can obtain HT on frequency- domain . 
The HT result can be deduced from FT, So 

x (f) = x(t)*h(t) = F~\x(f)H(f)] (3) 

As input signal x(t) = $m(2xf.t) , its FT 

X(f) = j 1 2[S(f + /.) - S(f - /. )] , and multiply the negative frequency of the X(f) by j and multiply positive 
frequency of the X(f) by -j 

X h (/) = \[Slf + f t )- S(f - ./,.)] = F[-co S (2^././)] (4) 

Take Fourier inversion to both sides of the formula 4: 

ff[siii(2jr/.f)] = -cos(2jr^.f) (5) 

Formula (5) shows HT of the sine function is the negative cosine function, the result also shows that it makes 
-j phase- shifting to the input signal. The same argument if input signal is the cosine function, the HT of the 
function is the sine function, the result shows that it makes -j phase- shifting to the input signal too. So the 

HT can provide 90* phase-shifting and not influence the amplitude of the frequency-spectrum components 
[3]. 

III. Basic Principle Of Phase Shifting Detection 

The ideal signal of power system is the pure sine waveform 

u(t) = Asin(2w/Qr + (Pq) = AsinOyr + (p^) (6) 

In formula (6): A is the voltage amplitude; / is the basic frequency: 50Hz; q> is the initial phase. 
In ideal conditions, the properties of u(t) are 

2 2 
u (r) + u (t-TIA) 

2 2 
= [A sin(e> r + q> Q )] + [A sin(e> (f - T 1 4) + p Q )] 

= A 2 

(7) 
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When the detection output y(t): 

y(t) = u 2 (t) + u 2 (t-T IA)-A 2 (8) 

Apparently, if input signal is the ideal waveform the detection output y(t) is zero. However when the 
input signal contains disturbances, the expression is 

«(f) = Asin(»Qf + <pQ) + e(j) (9) 

o 

After 90 phase-shifting to the u(t): 



The detection output y(t): 



H (0 = u(t-TIA) = Asin(io Q t + <p Q -90 ) + e(f - 90 ) (10) 



2 2 2 

y(t) = u (t) + u (t-TIA)-A 

= 2Asin(eW + q>n) * e(t) 

(ID 



+2 A cos(fi) I + q>Q ) * e (t - 90 ) 
+ e 2 (t) + e 2 (r-90°) 



In the formula 1 1, Asin(cw + q^) and Acos(cw + q>n) are known, so 

y(0 =£(«(»)) (12) 

So when there are disturbances in the input signal, the detection output signal, the detection output certainly 
cannot be equal to zero. Furthermore e(t) will have different value when the disturbances input are different 
and detection output will have different waveform characteristics, if it is corresponding to the time when 
disturbances take place, it can assure the detection's real-time characteristic^]. 

IV. Simulated Distrubance And results 

Various Power Quality disturbances like voltage swells, voltage sags, voltage fluctuation, harmonics, inter- 
harmonics and transient oscillation have been simulated and their real time detection was successful done. 

A. Ideal waveform generation 

An ideal waveform is generated in mat lab of frequency 50 Hz, amplitude 1 , sampling frequency 6400 Hz 
and number of samples 1280. When an ideal waveform with no distortion is given to Hilbert detector, the 
simulation shows a straight line, i.e., zero distortion. The ideal waveform is shown in "Fig. 1" and the zero 
distortion output is shown in "Fig. 2". 




Figure 1 . Ideal Waveform 
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Figure2. Zero Distortion Output waveform 
B. Voltage Swells and Voltage sags waveform 

u(t) = (l±a(Li(t-t l )-Li(t-t 2 )))sm((o Q t) (13) 

In formula (13), «(r) = j^'^ is tne unit ste P functions, a = 0.1 □ 0.9 Is the voltage swells (voltage dips) 

amplitude; 0.57" < t ~ < 30T is the voltage swells (voltage dips) durations; t| is the voltage swells (voltage 

dips) initial time, t 2 is the end time.[4]-[5] Here, a voltage swell signal is generated, using LabVIEW, whose 
amplitude is 1.9 times the normal amplitude during disturbance. Simulation results show that voltage swell 
occurs during the samples 320 to 700. When the distorted waveform is fed into the HT detector, a swell is 
seen during the disturbance. HT accurately detects the amplitude, duration and time of the voltage swell. 
"Fig. 3" shows the voltage swell in the generated wave and "Fig.4" shows the output of the voltage swell after 
phase shift. 



voltage swell Plct0 




Figure3. Voltage swells generation 




Time 



Figure4. Output of Voltage Swell after HT 
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A voltage sag signal of amplitude 0.2 times the normal amplitude during disturbance is generated and fed 
into the Hilbert detector. The simulation results show that sag occurs during the samples 320 to 705. "Fig. 5" 
shows the voltage sag generation and "Fig.6" demonstrates the sag detection. 




Figure5. Voltage sags generation 
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Figure6. Voltage sag Detection 



C. Voltage Fluctuation 

u(t) = (1 + a sin($f» ()) * sin(<5> (14) 

In formula (14), « = 0.10 0.2 is the fluctuation range; a = 0.1 □ 0.5 is the fluctuation frequency relative 
coefficient. A voltage fluctuation signal is generated using the equation (14), where the value of a=0.170 and 
P=0.206. The sampling info contains sampling frequency (6400 Hz) and number of samples 1280. Simulation 
results show a periodically fluctuating wave. The timed delay is taken to be O.lsecs. "Fig.7" shows the input 
data for voltage fluctuation, "Fig. 8" shows the Voltage fluctuation waveform and "Fig.9" shows the Output 
waveform for Voltage fluctuation. [7] 
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Figure7. Input data for Voltage Fluctuation 
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Figure8. Voltage Fluctuation Waveform 
VWo,m Graph S ^^^^^^B Plotl TO 
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Figure9. Output Waveform for Voltage Fluctuation 



D. Harmonics 



u(t) = sin(e?_r) + a s'm(3 G> t) + a s'm(5G>^t) + smiloy^t) (15) 

A distorted harmonic waveform is created by adding 3 ld , 5 lh , 7th and 9 lh harmonics to the ideal waveform. 
Results for detected output wave are given below in "Fig. 10" and "Fig. 11". [6] 




FigurelO.The Generated wave with harmonics 
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Figure 1 1. Output wave with harmonics 
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E. Frequency fluctuation 

In frequency fluctuation generation, during interval 450 to 700, signal of three times the normal frequency is 
added to the signal. Simulation results show greater fluctuations during that interval in "Fig. 12" and 
"Fig. 13". 

freq. fluctuation P>°*° SB9 | 




Figurel2.The input wave for frequency fluctuation 
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Figurel3.The output wave for frequency fluctuation 



F. Notch 

A notch denotes a disturbance that occurs for a very short duration. Hence, a signal is generated with 
disturbance at only one sample. Its simulation result shown below depicts a disturbance at sample 420 in 
"Fig. 14" and "Fig. 15". 




Figure 15. Output Waveform for notch 
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G. Momentary Interruption 

Momentary interruption means that a signal is interrupted during that time interval, i.e., signal is zero during 
that time interval. After executing the code, simulation shows error between samples 350 and 480. Figure 16 
and Figure 17 depict the aspects of Momentary Interruption. [8] 

momentary interruption PlotO BBS 




Figure 16. Input Momentary Interruption Waveform 
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Figure 17. Output Momentary Interruption Waveform 



V. Conclusions 

The phase shifting property of Hilbert Transform has been employed to detect and classify various faults in 
this paper. Disturbances such as voltage swells, voltage sags, voltage fluctuation, harmonics, inter-harmonics, 
transient oscillation, frequency fluctuation, notch and momentary interruption have been simulated and 
results have been found to be accurate with respect to detection of the faults, finding out their characteristic 
magnitudes and time of their occurrence. The method is hence suitable for Real-time application for the 
detection of faults as the time of occurrence is precisely detected. 
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